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Abstract. The construction of a Laplacian on a class of fractals which includes 
the Sierpinski gasket (SG) has given rise to an intensive research on analysis on 
fractals. For instance, a complete theory of polynomials and power series on SG 
^SJ ■ has been developed by one of us and his coauthors. We build on this body of work 

to construct certain analogs of classical orthogonal polynomials (OP) on SG. In 
^ ■ particular, we investigate key properties of these OP on SG, including a three- 

term recursion formula and the asymptotics of the coefficients appearing in this 
recursion. Moreover, we develop numerical tools that allow us to graph a number 
of these OP. Finally, we use these numerical tools to investigate the structure of 
■ the zero and the nodal sets of these polynomials. 
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1. Introduction 

The polynomials 7ro(t) = 1, and 7rfe(t) = + Yl^Zo '^'^^ k = 1,2, . . . defined on the 
unit interval / are called monic orthogonal polynomials with respect to a weight w{t) 
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if 

|2; 



(1) (vTfc, -Kj) = J 'n:k{t)'Kj{t)w{t)dt = HvrfcHs^j A; = 0, 1, 2, ... , 

where Sk,j denotes the Kronecker 6 sequence. 

When renormahzed, these monic orthogonal polynomials give rise to an orthonor- 
mal system of polynomials {vrfcl^g, where {ifk,7Tj) = 6j^k 

The monic OP are constructed by performing the Gram-Schmidt process on the 
monomials {t'^j^O' with respect to the inner product (a, b) = Jj a{t)h{t)w{t)dt. Ex- 
amples of such polynomials include the Legendre polynomials defined on / = [—1, 1], 
and using the weight function w{t) = 1, 

Numerical information about the OP can be obtained via a fundamental recursive 
equation called the three-term recursion formula. For example, the values of the 
polynomials at different points in J, and their zeros can be obtained from the three- 
term recurrence relation. For comparison with our results, we state the general form 
of this recurrence relation and refer to [HI [I3] for details on OP. 

Theorem A. [6, Theorem 1.27] The orthogonal polynomials {7!'k}k>o satisfy a 
three-term recurrence relation. More specifically, vr_i(t) = 0,7ro(t) = 1, and for all 
k > we have: 

(2) 7rfc+i(t) = (t- ak)nk{t) - /3fe7rfc_i(t), 



where 
(3) 




itei for k > 
, for k > 1. 



Our goal in this paper is to construct and investigate the properties of the analogs 
of the above OP on SG. But first what is a polynomial on a fractal set such as SG7 

Let A = ^ denote the Laplacian on / = [0, 1]. If we let vrj(t) = P ,foi t E I = [0, 1], 
it is immediate that ttj is a solution to a differential equation A^TTj{t) = 0, for some 
i > 1. This observation can be used to define polynomial on fractals in general and 
on the Sierpinski gasket in particular. In particular, on SG, one defines a polynomial 
P to be any solution of A^^^P = for some j > 0, where A is the (fractal) Laplacian 
to be defined below. Using this analogy, monomials were introduced on the Sierpinski 
gasket [H [m [21 [13 US] based on the fractal Laplacian constructed by Kigami |9], 
see also [151 [HI- We also note that a theory of polynomials based on a different 
Laplacian, the so-called energy Laplacian, is developed in [T7] . 

Consequently, our construction of OP on SG is based on these polynomials. In 
particular, we apply the Gram-Schmidt orthogonalization algorithm to produce the 
analogs to the Legendre polynomials on SG. In the process we construct both a 
family of symmetric and antisymmetric OP on SG, and the concatenation of the two 
families yields a tight frame for a proper closed subspace of L'^{SG). The reason 
that our construction does not yield a dense set of OP in L'^{SG) is due to the fact 
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that the polynomials are not dense in this space [9l Theorem 4.3.6], [TS", Section 
5.1]. In addition, we investigate thoroughly a corresponding three-term recursion 
formula, paying a particular attention to the asymptotics of the coefficients involved 
in this recursion. Our three-term recurrence is fundamentally different from ([2]), 
and this is essentially due to the fact that the product of two polynomials on SG 
is not a polynomial [1]! In addition, our analog of ([2]) is very unstable and thus 
cannot be used to generate the values of the OP on SG. However, we are able to 
exploit this three-term recurrence to recursively express the corresponding OP in 
terms of a well-known basis of polynomials on SG. This leads to the development of 
some numerical tools that allow us to plot graphs of many of these OP. In addition, 
we present preliminary results pertaining the zeros of these OP. In particular, we 
graph the zero sets of certain of these OP. These results are mainly experimental 
but seems to indicate that these zero sets have some structure. We have not been 
able to prove anything along these lines but hope that our experimental results lead 
to more investigations on the zero sets of the OP. Nonetheless, we have generated 
many graphs related to the OP on SG that we did not include in the present paper 
due to space constraint. The interested reader is referred to the following web site 
which contains many figures and algorithms that were generated in the course of this 
research project |19j. 




Our paper is organized as follows: In Section [2] we review the basic facts of analysis 
on fractals needed to state and prove our results. We also prove some new results 
about the Green's function that may be of independent interest. Section [3] is devoted 
to the construction of the OP on SG and to the investigation of the three-term 
recursion and the asymptotic analysis of its coefficients. We also consider the Jacobi 
matrix associated with these coefficients. In Section H] we carry out the numerical 
computations that enable us to plot many of the OP we constructed. In addition, we 
give some numerical evidence concerning the asymptotics of the coefficients involved 
in the three-term relation. Finally, we give numerical description of the zero sets of 
the OP and display some of these zero sets as well as some of their corresponding 
nodal sets on SG. 



2.1. Preliminaries. For more information on the theory of calculus on fractals see 
[9j; more specifically for calculus on the Sierpinski Gasket see [151 [H]. However, we 
collect below some key facts needed in the formulation of our results. 

The Sierpinski Gasket (SG), shown in Figure [H is the attractor of iterated function 
system (IFS) consisting of three contractions in the plane : — )■ M^, 2 = 0,1, 2, 
defined by = j{x — qi) + q-i where {gjj^^Q are the vertices of an equilateral triangle. 
SG satisfies the self-similar identity SG = L}f^QFi{SG), and the Fi{SG) are called cells 
of level 1. By iteration we can write = F^^^oFw^o- ■ -^F^^ for w = {wi, W2, . . . , Wm), 
each Wj = 0, 1, 2. We call each w a word of length \w\ = m, and we have the self- 
similar identity SG = U\w\=rnF'w{SG). Each F^^SG) is called a cell of level m. 
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Figure 1. The Sierpinski Gasket 

SG can also be viewed as the hmit of a sequence of graphs (with vertices Vm 
and edge relation x ~m y) defined inductively as follows: Fq is the complete graph 
on Vq = {qo, qi, q'2}, and Vm = Uf^QFiVm-i with x U if and only if x and y belong 
to the same cell of level m. Then K = ^m=i^ra-i the set of all vertices is the analog 
of the dyadic points in [0, 1] and is dense in SG. We consider Vq the set of boundary 
points of SG, and K \ Vq is the set of junction points. Note that every junction point 
in Vm has four neighbors in the graph F^. 

The graph Laplacian is defined by 

(4) {u{y) - u{x)) for a; e V;„ \ V^, 

and the Laplacian A on SG is defined as the renormalized limit 

(5) Au{x) = lim -b'^Amuix). 

m— >-oo 2 

At the three boundary points of SG, we have two different derivatives, the normal 
and tangential derivative, respectively defined by 

dnu{q,) = lim^^oo(i)"(2tx(g,)-ii(i^/"T/,+i)-«(^"*?i-i)) 
dTu{qi) = limm-,^5^{u{F^qi+i)-u{F^qi^i)), 



(6) 



identifying the boundary points Qi and qj when i — j mod 3. 

Using the above definitions, the analog of the Gauss-Green formula takes the fol- 
lowing form: 

(7) / {uAv - vAu)dii = y2{u{qi)dnv{qi) - v{qi)dnu{qi)), 

where ji is the natural probabihty measure which assigns weight 3~"* to each cell of 
order m. 
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These analytical tools can now be used to solve the boundary value problem 

-Au = /, u\vo = 0, 

whose solution is given by 



u= G{x,y)f{y) d^i{y). 

JSG 

G is called the Green's function, and is given by 



G{x,y) = hm GM{x,y), where Gm(x,i/) = ^ ^ 

m=0 z'&Vm+l\Vm 

where g{z, z') is zero when z and z' are not in the same cell of level m + 1, 

9 /SA"" 

g{z, z) = — {- ] for 2: e Vm+l \ Vm 



50 V5, 

and q(z, z') = — ( - 1 for 2; ^ z' , z and z' in the same cell of level m + 1, 
^ ^ 50 V5/ 

and '^/'i'"'''^'' is the piecewise harmonic function satisfying ilji'^'^^\x) = 6zx and vanish- 
ing outside the (m + l)th cell to which z belongs. 

We shall need an exact estimate of ||G||2,2, the LP' norm of the Green function G. 
It was conjectured in [10] and proved in [18] that 

= 178839/902500 ~ .198. 

This immediately implies that ||G||l2 < HGH/^oo < 1. However, we obtain exact values 
for llGllia and a related quantity in the next result. Though the result seems simple, 
to our knowledge, it has not appeared anywhere in the literature. 

First, we briefly recall a description of the spectrum of A on SG, that was given 
in [5] using the method of spectral decimation introduced in [12]. In essence, the 
spectral decimation method completely determines the eigenvalues and the eigenfunc- 
tions of A on SG from the eigenvalues and eigenfunctions of the graph Laplacians 
Am. More specifically, for every Dirichlet of —A on SG, that is a solution of 

—Am = XjU, M|yQ = 0, 

there exists an integer m > 1, called the generation of birth, such that if m is a Xj- 

(k) 

eigenfunction and k > m then u\v^. is an eigenfunction of A^ with eigenvalue A^- . 

The only possible initial values A^™''' are 2, 5 and 6, and subsequent values can be 
obtained from 



(8) Af = Wjz< fo, > ^ 

where can take the values ±1. The sequence A^'^^ is related to Aj by 

(9) A,- = I lim 5''Xf\ 



and 
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Theorem 2.1. Let {Aj}°^Q be the eigenvalues of the Laplacian on SG. Then, 
\\G\\h = / / G{x,yrdxdy = J2^ = Jflo ^ -0127 

JJSGxSG 

» CO 

/ ^(x,^)^^ = Vf = 1/6 

JSG ' 

Proof. It is easily seen that G can be expressed as 

oo 
k=0 

where are the Laplacian eigenvalues corresponding to the eigenfunction Uk [T5] . 
It then follows that 



[ G(x, x)dx = y f = lim I y" ^ 



(m) ; 
3 



where is a finite set that depends on the multiplicity of the graph-Laplacian 
eigenvalues A^-™'''. We shall now find explicit formula for the above sum. 

Suppose that we know the graph-Laplacian eigenvalues {A^'"'*} whose multiplicity 
are {d^^^}: if A^*"^ = 6 it has multiplicity if A^™^ = 5 it has multiphcity 

and if A^™'' = 3 it has multiplicity . Then, at the next level (m + 1), the 

graph-Laplacian eigenvalues {A^™'^^''} come from the spectral decimation method. 
Moreover, for every A^'"'' 7^ 6 with multiplicity d^™'' there will be two eigenvalues 
jj^(m+i) ^-^^Yi multiplicity d^J^^ and given by 



(m+l) _ 5±^25-4A^ 



(m) 



Note that 



^{m + l)_|_^{m + l) 



_l^{m + l) "T g„ + i_)^(m + l) 5m 5 ^(m + l)^(m + l) 



. (m) 



Now = li^^m^oo where Am = Ea,^6 ^^^^(m) . since the eigenvalue Xf"^ = 6 

with multiplicity contributes to the sum. Therefore, by the spectral 

decimation, we have 

nn^ A . - ^ s^+s 1 I 2 3'"'i-3 1 I 4 

V-"-"/ ^m+1 3 2 S^+^S 3 2 5'"+-'-3 

with = ||(l/2 + 1/5 + 1/5) = 3/25. Consequently, 

00 00 



A 



ORTHOGONAL POLYNOMIALS ON THE SIERPINSKI GASKET 7 

Similarly, if we let 5^ = | EA^6(^;;^)^' ^l^en 

= lim Brn- 

j 

One can check that 

{m+l),2 , ^, (m + l)^2 r,rr o\('") 

^™+lA*.™+^' 5m+i;^(™+i) 25 52'" (;^(™+l);)^('"+i))2 25 (5m;^{'"))2 25 ' 52m;^('") ' 

Next, we can use the spectral decimation method to prove that 

(^^\ R _ 4 3'"+3 1 I 4 3"'-3 1 I R _ Ar" 

I,-'--'-/ -Dm+1 g 2 ^5m+i .5)2 ~r 9 2 (5m+i.3'j2 ~r -Dm 25 



where = Ea,^6 s2J;,('") = §5^^™- Moreover, 



r — 3'"+3 1 I 3'"-3 1 I Ir^ 

'-'m+1 — 2 52m + 2.5 ~r 2 52m + 2.3 T g'-^m; 

with Ci = ^(1/2 + 1/5 + 1/5) = 9/250, and = §^(1/4 + 1/25 + 1/25) = 
Consequently, 



oo oo 



El — R _|_\^ 2 3"'-^+3 I 2 3'"-i-3 2 \ ^ ^ _ _l_ / 3532 167 \ _ 45389 ^ 97 

Af -°1"^Z^ 3 53-52™ "^Z^ 3 152-52'" 25 Z^^™ 25^52-81-11 1760/ 3564000 " '^"^^ ' ' 

□ 



j m=l m=l m=l 



2.2. Bases of polynomials on SG. For any integer j > 0, the set of polynomials 
of degree less than or equal to j will be denoted W and consists of the solutions of 
A^~^^u = 0. W is a space of dimension 3j + 3 and it has a basis {fki, < k < j;i = 

0. 1.2} characterized by 

(13) A'^fkiiQi') = 5e,kSi,i', 

where i' = 0, 1, 2 and < £ < j. In particular, "H" is the space of harmonic functions, 

1. e., a function h : SG — )• M belongs to "H" if and only if it satisfies Ah = 0. We refer 
to [m [m [15] for details on polynomials on SG. 

Our construction of OP on SG will be based on another basis for the space W 
introduced in [11]. Functions in this basis are called monomials and are essentially 
the fractal analogs of jf. 

Definition 2.1. Fix a boundary point g„ for n = 0,1,2, The monomials P^-"^ for 
i = 1,2,3 and j > are defined to be the functions in W satisfying 



(14) <^ d^A^P}f>{qn) = 

[ dTA-^P^\qn) = 5mj5,3, 
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with < m < j. When n = we will sometimes delete the upper index and just 
write Pji for Pji^ ■ 

Notice that for a fixed n = 0, 1, 2, the set of monomials {Pki, "^ = 1,2, 3};[^g forms 
a basis for W . Moreover, the monomials in this basis satisfy APj"^ = P^^J^^^-, and 

possess some symmetry properties. When i = 1 or 2, Pjf) is symmetric with respect 
to the line passing through g„ and the midpoint of the side opposing g„. In fact, these 
symmetric polynomials can be viewed as analogs of the even polynomials t"^^ on the 
unit interval. When z = 3, Pj^^ is antisymmetric with respect to the line passing 
through Qn and the midpoint of the side opposing g„. In this case, Pj^^ can be viewed 
as analogs of the odd polynomials t"^^^^. 

In our construction of OP on SG, we will need the following result which was 
proved in [TT] and which gives recursively values of the above monomials and their 
derivatives at boundary points. 

Theorem B [11, Theorem 2.3] For j > 0, let 

Oij = Pjiiqi), Pj = Pj2(gi), 7i = Pjsiqi), Vj = dnPjiiqi) 



The following recursion relations hold: 
(15) 



^ "i = 5^YllJiOij-eai for j > 2 



= T5(^ES(3-5^--^-5^+i + 6K_,/3, for j > 1 



5J-5 

15(5J 

7i = 3aj+i for j > 1 

. Vj = ^Oij + '^mZlvef3j-e for j > I, 

where the initial values are: = l,ai = 1/6, (Sq = —1/2, tjq = 0, and (9„Po2(Q'i) = 
—1/2. Moreover, 

dnPj2{qi) = -oij for i>l and dnPjsiqi) = ^Vj+i for j > 0. 

We can use Theorem B along with the Green-Gauss formula to compute the inner 
product among the monomials P^"'*- 

Lemma 2.1. Consider the monomials {Pji}j>o fori = 1,2,3. Then the following 
formulas hold: 

{Pji, Pki) = 2 Yle=j-m^ o;j-i7]k+e+i — c^k+e+iVj-e 
{Pj2, Pk2) = -2 Y.i=j-m, l^j-eOik+e+1 - l3k+i+iO!j-e 
{Pj3, Pks) = 18 Y,e=j-m, Oij-e+iVk+e+2 - ak+£+2Vj-£+i 
{Pji, Pk2) = -2 I^Lo '^'j^e'^'k+e+i - h+i.+2rij-i+i, 

where a'^. = at for k ^ 1 and a[ = —1/2, and where m^, = min (j, k). 

Moreover, 

(17) {PjuPk-i) = {Pj2,Pk3)=0. 



(16) 
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In addition, for i = 3 the inner products among the different monomials associated 
with each of the boundary points g„, n = 0, 1, 2 are given by 



[18) 



for n ^ n' . 



/p(") p(")\ _ /p(0) p(0)\ 

\-^j3 ' -^fcS / ~ \-^i3 ' ^k3 I 

/pin) pK)\ _ _i/p(0) p(0)\ 

\-^i3 5 ^k3 I ~ 2 \-^i3 ' ^kZ I 



Proof. It is enough to prove that for each k,j and any we have: 



i 2 

(19) (Pji, Pki') = ^ ^ {P{j-e)iiqn)dnP{k+l+e)i'{qn) — P{k+l+£)i'{<ln)dnP(j~e)i{qn)) 
£=0 n=0 



Notice that by a symmetry argument we can assume that k > j and proved the result 
by induction on j. Fix i,i' G {1,2,3} and notice that for all A; > Green-Gauss' 
formula gives: 



{Poi.Pki') — / PoiPki' dfi — / PoiAP(^k+i)i' dfi 
JsG JsG 

2 

APoiPki' dfl + Poiiqn)dnPki'{Qn) — Pki' {Qn)dnPoi{qn 

n=0 

Poi{qn)dnP(k+l)i'{qn) - P{k+l)i' {qn)dnPoi{qn) , 



n=0 



where we have used the fact that APoi = 0. This establishes ( !T9|) for j = 0, and all 
A; > . 

Now assume that for some j and for k > j we have 



j 2 

{Pji, Pki') = {P(j-£)i{qn)dnP{k+l+£)i'iqn) — P{k+l+t)i'{qn)dnP{j~t}i{qn)) , 

£=0 n=0 
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we just have to show that this last relation holds for j + 1 and k > j + 1. But, by 
the Gauss-Green formula we have 

{P(j+i)i, Pki') = / P{j+i)iPki' = I P(j+i)iAP(^k+i)i' dfi 

JSG JSG 
2 

= {Pji, P{k+l)i') + ^ {P{j+l)i{qn)dnP{k+l)i'{<in) " P{k+l)i'{<ln)dnP{j+l)i{<in 
n=0 

i 2 

= ^ ^ {P{j-t)i{qn)dnP{k+2+t)i'{qn) — P(k+2+t)i'{qn)dnP{j-t)i{qn)) 
l=Q n=0 
2 

+ ^ {P(j+l)Mn)dnP{k+l)i'{qn) — P{k+l)i' {qn)dnP{j+l)Mn)) 
n=0 
i+1 2 

= {P{j+l-t)i{qn)dnP{k+l+t)ii{qn) — P{k+l+e}i>{qn)dnP{j+l-t)i{qn)) 

e=i n=0 

2 

+ ^ {P(j+l)i{(ln)dnP{k+l)i'{qn) " (5'n)c^n-P{j+l)j(q'n)) 
n=0 
i+1 2 

= (-P(j+l-^)«(?ra)'9n-P(fc+l+^)i'(Q'n) " P{k+l+£)i' {qn)dnP{j+l-£)i{qn)) 

e=o n=0 

This shows that ( !T9|) holds for all j > and k > j. 

Next, f|T6|) follows from f|T9|) by symmetry consideration and by observing that the 
values and normal derivatives at go vanish. In addition. Theorem B, especially ( fTSl) 
can be used to further simplify ( |T9l) and establish ( fT6|) . 

Note that ( fTTl) is trivial because of the symmetry of Pji and Pj2 and the anti- 
symmetry of Pj3. 

To prove f fTSj) . it is enough to show that 

/p(0) p(l)\ _ _l/p(0) p(0)\ 

since the other inner products can be computed using similar arguments. 
Using Equation (fT9l) . and canceling out zero terms yields 

3 

{Pjl ^Pkl) = X] ^(°-£)3(52)5n^'(fc+£+i)3(g2) - ^(fc+£+l)3 ('?2)<9n^'(j%3 (52) 
£=0 

j 
^=0 



_l/p(0) p{0)v 



□ 
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3. Orthogonal polynomials on SG 

We are now ready to construct the families of orthogonal polynomials on SG. These 
OP will be obtained by applying the Gram-Schmidt algorithm on the polynomials 
{Pjk^}j>o, where k = 1,2,3, and n = 0, 1, or 2. In Subsections 14.11 and 14.21 we plot 
the graphs of various OP obtained from {Pjk}j>o, and k = 1,2,3. As mentioned in 
the introduction, more graphs of these OP can be found at [T9] . 

But we first derive a three-term recursion formula for our OP and give estimates 
on the size of the coefficients appearing in this recursion formula. These coefficients 
are extremely small. This creates major instability problems when plotting the corre- 
sponding OP. To circumvent this problem, we carry out our numerical simulation to 
arbitrary precision. A similar phenomenon was already observed in [11], and rational 
arithmetic was used. We also considered the use of rational arithmetic, but because 
it involves a prohibitive time cost, we settled for arbitrary precision instead. 

3.1. General theory of orthogonal polynomials on SG. 

Definition 3.1. Fix A; = 1, 2 or 3 and denote by {pj}^o {Pjk}'jLo the orthogonal 
polynomials obtained from {Pjfc}°^Q by the Gram-Schmidt process, i.e., pj = Pjk — 
^iZo d'i{Pjk,Pi)pi for each j > 1. Consequently, there exists a set of coefficients 
{uj^i}j^Q, with LOjj = 1, and such that 

j-i i-i 

(20) pj{x) = Pjk - ^ df{Pjk,Pi)pi = Pjk{x) + ^ ujj^iPikix),j > 1. 

/=0 1=0 

Moreover, 

{Pj,Pi) = dj^Sj^iwheie \\pj\\l2 = dj^. 

By normalizing the orthogonal polynomials {pj}'jLo obtain the family of or- 
thonormal polynomials {Qj}^Q characterized by 

{{Qj,Qk) = 5j,k 
Qj = djPj = djPjk + dj YjIZI ^j,iPik,j > 1- 

Theorem 3.1. Given k = 1,2 or 3, and for each j > the following holds: 

WpjUl^ = df^ < ||Pj,fc||L2. 

Moreover, for any < r < oo, there exist constants Ci, c,- > such that for all j > 

(21) WpjWl^ = dj' < Ci(j!)"'°s5/log2 ^ 

In particular, 

lim IIPiHl^ = lim dj^ = 0. 
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Proof. Recall that by the definition of pj we can write that for j > 1, the Gram- 
Schmidt process yields pj = Pjk — ^iZo d'i{Pjk,Pi)pi for each j > 1. Therefore, 

i-i i-i 

1=0 1=0 

which establishes the first estimate. 

Now recall from [HI Theorem 2.7] and [TTl Theorem 2.13], that for each r G (0, oo), 
there are c > and > such that 

0<aj< c(j!)-i°g5/iog2 Q < |^^.| < ^^^-j J > Q_ 

When k = 1, we get from Lemma [2. II that 

j 

\\Pji\\h = 2 "i-'^i+'+i - (^j+i+iVj-i 

1=0 

j j 



2 1 ^(yir]2j^i+i - y^a2i-;+i^;| 



j j 



< 2cCrr-'-^ ^(;,)-log5/log2 ^ 2cC,((j + 1) l)" 1°S 5/ log 2 ^ 



log 5 / log 2 

i=0 «=0 



from which fl?I]) follows. When = 2 or = 3 similar arguments are used to obtain 
the same estimate. It then follows that limj_j.oo dj^ = with at least an exponential 
rate of decay. □ 

Our first main result is the following theorem that establishes the three-term re- 
cursion formula for OP on SG. 

Theorem 3.2. Let {pfcj^Q be the orthogonal polynomials defined above. Let fo{x) = 
and for k >0, let fk+i be the polynomial defined by 



(22) fk+i{x) := - J G{x,y)pk{y)dij{y). 

Set p^i{x) := 0, and po{x) = Pq^^x). Then, for each k>0 

(23) Pk+iix) = fk+iix) - hkPk{x) - CkPk-iix), 
where 



dk{fk+l,Pk), 

Consequently, 



(24) ^ _ .Xfc.i _ iiF/cii^a 



(25) d^. = WpkWii = do C1C2C3 . . . Ck-iCk 
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Proof. Using the definition of fk+i we can write 

{fk+i,Pj) = - G{x,y)pkiy)pjix)dfx{y)dfx{x) = {fj+i,Pk) 



Notice tliat /j+i is a polynomial of degree j + 1. Tfius, /j+i is orthogonal to all pk 
with j + 1 < k. Therefore, fk+i = cLkPk+i + bkPk + CkPk-i for some coefficients a^, bk, 
and Cfc. It is easy to see that = 1 this follows from the fact fk+i —Pk+i is orthogonal 
to pk+i- Thus, fk+i = Pk+i + bkPk + CkPk-i- Taking inner product with pk yields the 
first equation in ^IM . Taking inner product with Pk-i yields 



Ck{pk-i,Pk-i) = Ckd^\ = {fk+i,Pk-i) 

Gix, y)pk{y)Pk-i{x) dnix) d^{y) 



'SG JSG 

= {fk,Pk) = {Pk + bk-lPk-1 + Ck-lPk-2,Pk) 

= df 

which is exactly the second equation in (^^. 

Observe that the last equation in is equivalent to Ckd^ = d\_^, which im- 
plies ([25]). □ 



The above results deserve some discussions. While fl23|) resembles ([2]), these two 
relations are fundamentally different. In fact, ([2]) is essentially the statement that 
tTTfc is a monomial of degree /c + 1, and thus can be expressed in terms of tt;, / < + 1. 
Because the product of polynomials is not a polynomial on SG, tn^ is replaced by 
fk+i{x) = — J G{x, y)pk{y) dfj,{y). Therefore, unlike in the classical case in which all 
information on Tik+i can be gathered from the coefficients in ([2]), on SG we must 
also evaluate the auxiliary polynomial fk+i- This presents an additional difficulty in 
carrying out any numerical simulations with our OP. 

We now prove a version of Theorem 13.21 dealing with the three-term recurrence for 
the orthonormal polynomial {Qk} 



oo 

fc=0- 



Theorem 3.3. Let {Qk}kLo be the orthogonal polynomials defined above. Let /q 
and for k > 



(26) fk+iix) :=- j Gix,y)Qkiy)dfiiy). 
Then, for each k >0 

(27) y/Ck+iQk+i{x) = fk+i{x) - bkQk{x) - y/(^Qk~i{x), 

where Q-i{x) := 0, Qo{x) = d^Posix), and b^. and Ck were defined in Theorem \3.2[ 

Proof. The proof follows from the fact that pk = \\Pk\\L'^Qk and (??). □ 
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We prove below certain properties about the coefficients bk and appearing in the 
three-term recursion formula. 

Theorem 3.4. For each k >0, we have 

(28) - \\G\\l2 <bk<0, 

(29) < Cfe < 
and 

(30) d,l, < \\Gh.d,\ 
In particular, 

< c?o^II^IIl2' 

1/2 



where \\G\\l2 = ^ ffgQ^g^G{x,yYdxdy^ 



Proof. \bk\ = dl\{fk+i,Pk)\ < dl\\fk+i\\L4Pk\\L^ < dl\\G\\L4Pk\\l2 = \\G\\l2. 

Using the expression of G{x,y) given in the proof of Theorem I2.H we can write 

oo 

(fk+uPk) = - ^X^^{pk,ui)^ < 0. 
1=1 

Hence, bk < 0. Similarly, 

< Ckd^\ = {fk+l,Pk-l) < ||/fe+l||L2|bfe-l||L2 < ||G'||L2||pfc||i2||]9fe_l||i2 = ||G||i2rf^M^\. 

Using the fact that Ck = d^'^dl^^, we have 

< Cfc < \\G\\L2dk-id'^'^ = \\G\\L2y/(^, 

from which the second estimate follows. 
The last estimate follows by observing that 

< d,l < Wfk+iWh = dll, + bldf + cldl\ < WGWUf. 
Consequently, 

dl'<d-,'\\G\\l2. 

□ 

Remark 3.1. The last estimate in Theorem 13.41 along with the estimate on the ||G||l2 
given in Theorem 12. gives another proof that the sequence d'^^ decays exponentially 
fast. 
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3.2. Recurrence relations for orthogonal polynomials on SG. We have now 
proved via two different arguments that the coefficients d^^ are very small. This, 
together with the fact that we must also evaluate the auxiliary polynomial fk+i 
prevents us to use Theorem 13.21 to directly plot the OP we construct. Rather we 
use (123|) together with the representation of pk in Definition 13.11 to plot these OP. 
More specifically, we identify each pj with a vector Qj G W"^^ where = ujj^i for 
/ = 0, 1, . . . , j — 1 and = 1. This can be simply written as 



( ^j,o \ 



( ^j,o \ 

\ 1 J 



We now derive a recursion relation for the coefficients fij appearing in the above 
representation of pj. We shall later, implement this recurrence using arbitrary preci- 
sion to give accurate graphs of Qj = dj^pj. 

Theorem 3.5. The recurrence relation for the orthogonal polynomials pk is accom- 
plished by evaluating the following equations for each k > 0. 



or more specifically, 



(31) 



^k+l,k 
1 



Ck 

^k,0 



^k,k-2 
^k.k-l 
1 



^k,0 
^k,l 

^k.k-l 
1 





- Ck 



( ^k-l,0 \ 
^k-l,l 

^k-l,k-2 
1 






where bk, Ck where defined in Theorem \3.2\ and (k = '^e=o^k,£'~fi+i, with 'ji given 
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Proof. Using fl22|) we can write 

j 



1 



G{x,y)pj{x)dn{y) = - / G{x,y)^ujj^iPe^kdn{y) 



i=0 



£=0 \ J / i=o V 



£=0 

j 

f \ 

\ 1 / 



7m p 
70 



□ 



We now derive a version of the Christoffel-Darboux formulas for the OP {Qk}- 
As an apphcation, we use these formulas to find the expansion of AQk in terms of 

{Qi}i=o- 



Theorem 3.6. Let N > 0, x,y E SG, and define K]\f{x,y) = Yl!k=oQk{x)Qk{y)- 
Using the notations of the last theorem we have 



N 



(32) K^ix^y) = y/cN+i[QNix)AQN+iiy) -QN+iix)AQNiy)] + '^ fk+i{x)AQkiy). 
In particular, for each x G SG and N > 0, 

(33) Kn{x,x) = y/CN+i [Qn (x) AQat+i (x) - Qn+i (x) AQjv ix)] + '^ fk+i {x)AQk {x) . 



k=0 



N 



k=0 



Proof Let A; G {0, 1, ... , N}, and x,y e SG. By (ETj) we have 

fk+i{x)Qk{y) - fk+i{y)Qk{x) = y/ci^[Qk+i{x)Qkiy) - Qk+i{y)Qk{x)]- 

Vci^[Qkix)Qk-iiy) - Qk{y)Qk-iix)]. 

Therefore, summing both sides we obtain 

N 

(34) ^[fk+iix)Qkiy) - fk+iiy)Qkix)] = y/cN+i[QN+i{x)QN{y) - QN+iiy)QNix)]. 

k=0 
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Now, observe that 

N N ^ N 

^[fk+iix)Qkiy) - fk+iiy)Qkix)] = ^fk+iix)Qk{y)+ / ^G{y, z)Qk{z)Qk{x)dz 

k=0 k=0 "^^^ k=0 

N 



^fk+i{x)Qk{y)+ G{y,z)KN{z,x)d2 

I — n J SG 



k=0 

Using this last equation and taking the Laplacian with respect to y of both sides 
of ([31]) yields, 

^lfk+ii.^)Qk{y) - fk+i{y)Qk{x)] j = ^ycN+l [Qn+i {x)AQn (y) - QN(,x)AQN+i{y)) 

^ k=0 ^ 

N 

= ^ fk+iix)AQkiy) - Kn{x, y) 

k=0 

Consequently, 

N 

KN{x,y) = y/cN+i[QN{x)AQN+iiy) - QN+i{x)AQN{y)] + fk+l{x)AQk{y)■ 
k=o 

□ 

We now use this theorem to find the coordinates of AQk in terms of lower order 
polynomials. That is since, AQk G span{Qi}'j^~Q, we have 

fc-i 

AQk{x) = Y,A^i''^Qii^), 

1=0 

where a[^^ = {AQk,Qi). These coefficients can be computed recursively using The- 
orem [321 More specifically, we have 

Corollary 3.1. Let k > 1, and assume AQk{x) = Yle=o ^^''^Q^i-''^)- Then, A^^}-^ = 
and for £ = 0,1, ■ ■ ■ ,k — 2 we have 



Ck ■ 

(fc) _ 



with the initial condition AQi = A^q^'Qq where aI^^ = -j=. 



Cl 



Furthermore, for each k > 1 and i = 0,1, . . . ,k — 1, A^^^ > 0. 

Proof. Use = in ( !33l) . gives Ko{x,x) = Ql{x) = ^yclQo{x)AQl{x), and inte- 
grating this equality leads to 1 = y/ci{AQi,QQ) . Hence, AQi{x) = A^^Qo{x) = 

:^Qoix). 

Next, write AQ2{x) = A^^^Qi{x) + A''^^Qo{x). Using = 1 in ([33D gives 
Ki{x,x) = y/^[Qi{x)AQ2{x) - Q2{x)AQi{x)] + f2{x)AQi{x) 
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where we use the fact that AQq = 0. Integrating this last expression leads to 

2 = y^[(AQ2,Qi) - ^{Q2,Qo)] + ^{f2,Qo) = v^4'^ + ^V^i. 

Therefore, = ^=. 

To compute Aq we utilize ( 132|) with = 1. In particular, we multiply Ki(x,y) 
by <5i(a;)Qo(y) and integrate the resulting equation with respect to both x and y. 
This will give 



(2) _L J?^ 
Cl 



= y3^(AQ2,Qo) + ^{f2,Qi){Qo,Qo) = 7^4'^ + 

Consequently, A^^^ = — = -A^^^^. Observe that Mx) = ^Q2ix) + 
+ y^Qo(2;) was used in the above arguments. 

Assume now that we have determined the coefficients {^^'^''}^=o ^Qk with re- 
spect to the orthonormal set {Qe}eZo- We can use an induction argument to find 
the coordinates of AQfc+i in the orthonormal set {Qe}i=Q- Set N = k in (l5H]) . and 
integrate the resulting equation to obtain A^'^'^^^ — ^ 

To obtain the remaining coefficients, A^'''^^^ for / = 0, 1, ■ ■ ■ , /c — 1, set N = k 
in f l32|) . multiply the resulting equation by Qk{x)Qi{y) and integrate with respect to 
both X and y. 

To obtain the positivity of the coefficients A^^\ we proceed by induction on k. 
Clearly, A^^^ = -i= > 0, and = -A'^^^ > since bk < for all k. In addition, 

(2) 1 

= ^= > 0. For assume, that for k > 2, and each £ = 0, 1, . . . , /c — 1, we 

i 1 — 1 1111 

have aI'^^ > 0. Then, aI.^'''^^ = , ^ > and for each £ = 0, 1, . . . , A; — 1 we have 

A^*^^"^'' = —Af^^ y/ck+i ^ ^ since 6^ < and Af^ > by the induction hypothesis. This 
concludes the proof. □ 

Remark 3.2. Note that Theorem 13.61 and Corollary 13.11 hold for the symmetric OP 
Sk defined in Subsection 14.21 



We end this section with an investigation of the tri-diagonal symmetric Jacobi 
matrix corresponding to the OP {Qk}- As in the classical case of OP defined on 
the real line, each of the families of orthogonal polynomials we constructed above is 
related to a tri-diagonal, symmetric Jacobi matrix that we denote J„. In particular, 
we can write fl27p of Theorem 13.31 in a matrix form that will involve this Jacobi 
polynomial. To do this we use the following notations: 

= [foi 111 - ■ ■ 1 fn-lY' <5^"^ = [Qoi Qli ■ ■ ■ 1 Qn-l]'^i 
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and 



Jn 



( bo 







C2 
b2 














\ 






\ ■•• yc;rT &n-i / 

By writing ( 1271) in a matrix-vector form we have for each n > 2 

(35) = J„g(")(a^) + v^Qn(x)e„ 

where = [0, 0, . . . , 0, 1]^. 

For the OP we consider in this paper, though J„ is invertible for all n, its de- 
terminant is extremely small. Recall [SI Theorem 0.9.10] that the determinant of 
a tri-diagonal symmetric Jacobi matrix obeys the following recursion formula. Let 
Dn = Det{Jn), then 

Dn+l = bn-lDn — C„_i-D„_i 

with initial conditions Dq = b^ and D_i = 1. 

To illustrate the fact that the Jacobi matrices associated to the OP considered here 
have small determinant, we show below a graph of log as a function of n. This 
graph is generated using the antisymmetric OP that will be constructed explicitly in 
the next section. 



sequence of Del for Antisymmetric OP 




100 
k 



Figure 2. Graph of Log\Dk\ of the determinants of the Jacobi matrix, Dj, 



4. Numerical results 

We shall now generate graphs of the OP we constructed above starting from two 
families of monomials on SG. We also present graphs of the sequences appearing in 
the three-terms recursion formula. We start by considering a family of antisymmetric 
monomials. 
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4.1. Antisymmetric orthonormal polynomials. Recall the monomials {Pjs} in 
Definition 13.11 are centered around the point go and antisymmetric with respect to 
that point. We now use these monomials to generate the family of antisymmetric 
orthogonal polynomials {pj}j^Q and their corresponding normalized version {Qj}°^q. 
We do this inductively using Theorems 13.21 and 13. 5[ In the process, we compute the 
sequence dk,Ck and 6^. 

More specifically, it is easily seen that ||Po3||i2 = d^'^ = 3/10. Next we use Theo- 
rem 13.51 to find the auxiliary polynomial 



Notice that we impose the condition fj\vo = 0. Thus, bo can now be computed, 
which , together with the polynomial /i, is used to find the polynomial pi. Next, 
/2 is computed again using Theorem 13.51 from which one gets di. The next step 
is to compute Ci. Using this, bi and the polynomial p2 are calculated yielding C2. 
Continuing these recursions, we construct pk, k >2 and all the related sequences. 

To better display the results of these computations, we use a logarithmic scale to 
plot the numerical sequences obtained above. In particular, we plot log(||pfc||^2 ) = 
\og{dl), log(— and logc^ versus k. These graphs are displayed in Figure [31 In 
particular, these graphs illustrate Theorems 13.11 and 13.41 That is ||pfe||^| tends to 
infinity exponentially fast. In addition the graphs of \og{—bk) and logc^ suggest that 



for some ko. This seems to be the case for k < 200, but we do not have data beyond 
this range of k to confirm or infirm the observation. 

In Figure m we plot the orthonormal polynomials Qj = djpj for j = 0,3,4,7, on 
Fy the level 7 approximation to SG. To do this, we use the recursion formula given 
by Theorem 13.51 along with the second half of the representation of pj given in (120|) . 
as well as the graphs of the monomials Pj^3 which were given in [11] . Moreover, most 
of the computations have been carried out to arbitrary precision since the coefficients 
Uj^i are very small. We refer to Table [T] for approximate values of these coefficients 
for j = 0,1,...,6. 

The fact that Qk are antisymmetric is apparent from Figure HI In addition, we 
observe that for i = 1,2, \Qk{(li) \ is larger compare to for any x G SG\{qi, 52}- 

This is consistent with the behavior of classical OP such as the Legendre polynomial 
on [0, 1]; see [7]. 

4.2. Symmetric orthogonal polynomials. We next carried out the constructions 
of subsection 14.11 to the fully symmetric polynomials which we first define. By fully 
symmetric we mean symmetric under all rotations and flips in D^. Then we will 
apply the results of section [3] to this family of fully symmetric polynomials to obtain 
the symmetric OP {sytj^o their normalized counterpart {Sk}'^^^- 





Ck < —bk for k > k[ 



■0 
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Sequence of Log|rf,-| (Antisymmetric) Sequence of Log[-i;i|(Antisymmelric) 




200 



k 

Figure 3. Plots of the \ogdl, log(— 6^), and logc^ in log scale for the 
antisymmetric OP pk- 



Table 1. pk and 





Po,3 


^1,3 


-^2,3 


Po 


1 








Pi 


-2.04E-02 


1 





P2 


1.08E-04 


-1.30E-02 


1 


P3 


-3.23E-07 


6.03E-05 


-9.74E-03 


Pi 


1.99E-10 


-6.41E-08 


2.14E-05 


Pb 


-8.20E-14 


5.16E-11 


-2.94E-08 


P6 


4.43E-17 


-3.53E-14 


2.63E-11 



their coefficients ooji 



-^3,3 


^,3 


^5,3 


-^6,3 






































1 











-6.01E-03 


1 








1.41E-05 


-5.03E-03 


1 





-1.77E-08 


9.99E-06 


-4.23E-03 


1 



Definition 4.1. Define the fully symmetric monomial as follows for each j > 0: 
(36) P, = Pi?+Pi? + Pff 

The symmetric OP denoted {sj}^Q are now obtained from these symmetric polyno- 
mials by applying the Gram-Schmidt process and satisfies 

where dj'^ = Moreover, Sj = pj + J2iZo f^j,iPi for a set of coefficients {fij^iYilQ- 

By normalizing the orthogonal polynomials Sj we obtain the orthonormal OP de- 
noted {S-.j-o = {d,Sj}JL,. 



22 KASSO A. OKOUDJOU, ROBERT S. STRICHARTZ, AND ELIZABETH K. TULEY 

Remark 4.1. Note that the above definitions will remain unchanged if in defining the 
symmetric polynomials pj we used Pj2 instead of Pji; see for details about this. 

Observe that for each j,k > we have {pj, pk) = 6(Pj°\ -Pfc?); which implies that 
for j > 0, ||pj||^2 = dj'^ = 6||Pji||^2- We remark that the size of ||sj||i2 was already 
estimated in Theorem 13.11 Furthermore, one checks easily that d^^ = ||po|li2 = 
6||i'oi|li2 =6 

As in subsection 14. H we now plot the symmetric orthonormal polynomials as well as 
some of the related sequence. In particular, Figure [5l displays plots of the sequences 
log(||sfc||^2 ) = log((ifc), log(— and logc^ for these symmetric OP. The behavior of 
the coefficients Ck,bk, and d^ in the case of the symmetric OP Sk are very similar to 
those we observed for the antisymmetric OP. 

In Table [2], we show the coefficients pj^g for Sj, when j = 0, 1, . . . , 7 

Four of the symmetric orthogonal polynomials Sj = djSj, corresponding to j = 
0, 3, 5, and 6 are shown in Figures [61 As observed for the antisymmetric OP, \Sk{qi)\, 
which is constant for i = 1, 2, 3 (due to symmetry) is very large compare to the values 
of Sk at non-boundary points in SG. 

4.3. Orthonormal system of polynomials. By combining both the antisymmet- 
ric and symmetric OP, we can form an orthogonal system of polynomials in L'^{SG). 





Figure 4. 4 antisymmetric orthonormal polynomials. 
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Sequence of Log[rf,-j (Symmetric) Sequence of Log[-i>t| (Symmetric) 




k 

Figure 5. Plots of the logdl, log(— 6^), and logc^ for the symmetric 
OP Sk in log scale. 

Table 2. sj and their coefficients fij^e 





Po 


Pi 


P2 


P3 


P4 


P5 


P6 


So 


1 




















Si 


-5.56E-02 


1 

















S2 


5.34E-04 


-2.44E-02 


1 














S3 


-6.98E-07 


9.72E-05 


-1.23E-02 


1 











sa 


1.56E-09 


-3.20E-07 


6.15E-05 


-9.93E-03 


1 








S5 


-1.46E-12 


3.54E-10 


-9.41E-08 


2.62E-05 


-6.48E-03 


1 





S6 


1.44E-16 


-8.51E-14 


5.21E-11 


-2.94E-08 


1.41E-05 


-5.02E-03 


1 



Notice, that this system will not span the whole space L^, since the space of polyno- 
mials is not dense in L'^{SG) either. 

Using Lemma [2.11 we can prove that the set {Pj^\ Pjl\ P^"^^} forms a tight frame 
for its span J^j = span{Pj^\ Pjs^, Pj'^}- That is, there exists a constant Aj > such 
that for each P G J^j = span{Pj^\ Pjl\ Pj'^} we have 

2 

Pi^) = A,Y,{P,p\f)p\-\x),\/x e SG. 

n=0 

For more on frame theory we refer to [3]. 
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Theorem 4.1. For each j > 0, the set [pj^\Pjl\Pj^^^ forms a tight frame for its 



two-dimensional span Tj = span{Pj^\ Pj^', Pjs'} with frame bound A 



,(1) 



i3 



(0)| 



2 

L2- 



Proof. The frame bounds of [Pj^\ Pjl\ Pj^^ } are determined by the eigenvalues of 
its Gram matrix G given by 



G 



p(i) 



,(0) p(l) 



J3 



j3 



J3 , 



-IIPW||2 



2 -1 -1 
-1 2 -1 
1 -1 2 



G has eigenvalues {0, | 
operator 



|p(0)| 



L2' 2 



\Pj^^ Ilia}- This automatically shows that the frame 



,(0) 



j3 



j3 



□ (2) 

i3 



p(l) 

^i3 
yp(2) I 

V j3 / 



is the 2 X 2 diagonal matrix with fHP-J^Hia on the diagonal. Thus, {Pj^\ Pff, Pj^^} 
is a tight frame with frame bound Aj = MPj^^ 11^-2 □ 




Figure 6. 4 symmetric orthonormal polynomials. 
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Recall that {Sj}'jLQ are the symmetric orthonormal polynomials, and for each i — 

0,1,2, corresponds a family of antisymmetric orthonormal polynomials {Qj*''}^o- 
With these notations, the following result holds. 

Theorem 4.2. The set 



-Qr' + \-Sj,^ = 0,1,2 



3=0 



is an orthonormal basis for a subspace V of L'^{SG) . 

Proof. For j > and i e {0, 1, 2}, let an element in the above system be defined by 
Then, the inner product of two orthonormal system elements is: 



Recall that Sj is fully symmetric and Q^"'' is antisymmetric, so their inner product is 
zero for all j and n. Moreover, ii n ^ n',hy the Gauss Green formula, we can rewrite 

{Qt\Qt'^) -Y.Y.Qf\qi)dnQf\qi) - Qf\qi)dnQf\qi). 

j=0 i=0 

We now evaluate the right-hand side of the last equation. When i = n or i = 
n' , Q^^\qi) = 0, and this eliminates two terms for each j . Let Oj = Q^J^\qn+i), 

Then we have 



then —Qj = Qj (qn-i) (by symmetry). If we let bj = dnQj {qn+i), then —bj = 
dnQf'iqn-i)- H we let Cj = dnQf\qn), then = since dnP-^\qn) = for all j. 



k 
3=0 

More generally, the same argument can be used to prove that, {Q^^\ Q^? ^) = (^fe,it'(^n,n'- 
Consequently, 

□ 

Remark 4.2. Notice that using a different normalization we can show that the set 
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is also an orthogonal system for L'^{SG). 

Four of these orthonormal polynomials are plotted in Figure U\ and as might be 
expected, these OP do not seem to possess any obvious symmetry. 








Figure 7. 4 orthonormal polynomials. 

4.4. Zero sets of orthogonal polynomials on SG. Classical orthogonal polyno- 
mial theory suggests that zeros of higher order polynomials should interweave between 
the zeros of the next higher order polynomial. Specifically, if xq is a zero of pn, then 
there exist Xa and Xb zeros of P„+i such that Xa < xq < Xb- 

On SG, we have not been able to establish a similar result. In fact, it is not clear 
what 'interweaving' will mean in the fractal setting. More generally, on SG, the zero 
sets of the OP seem very difficult to fully describe. However, using ( 135|) we see that 
an exact characterization of the zero set of Q„ is the set of all points xq G SG such 
that 

Though this looks similar to the characterization of the zeros of classical orthogonal 
polynomials in terms of the eigenvectors of a Jacobi matrix, there is a main difference 
in that F„(xo) is a function of the auxiliary polynomials fk, k — 0,1, . . . ,n — l. Despite 
this difficulty, we have some numerical data, suggesting some structures in these zero 
sets. In particular, we shall display below nodal domains corresponding to some of 
these OP. 
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We first consider the zeros of tlie OP on tlie edges of SG. Notice tliat tliere are 
two types of edges, wfiicli we call bottom edge (the edge across from go) aiid side 
edges (edges which meet go)- We present some numerical data in figures [HI El andfTOl 
which seem to indicate that the zeros on these edges behave like zeros of classical 
orthogonal polynomials. Because 5*^ is fully symmetric its behavior on side edges is 
the same as their behavior on the bottom edge. The restriction of Qk to the side 
edge of SG as illustrated in Figure [TOl 





Figure 8. Graph Q^, k = 0,1,4,7 on the bottom edge of SG 




bottom edge 



bottom edge 



bottom edge 



bottom edge 



Figure 9. Graph of Sk, /c = 1, 4, 5, 6 on the bottom edge of SG 





side edge 




Q 



10 




side edge 



Figure 10. Graph of Qk, k= 1,5, 6, 10 on a side edge of SG 

More generally, we have been able to numerically generate the nodal domains of 
both the first 20 antisymmetric and first 20 symmetric OP on SG. For example, the 
nodal domains of a polynomial / on SG are defined as follows: Let 

Z{f) = r'{0} = {xeSG: fix) = 0} 

be the zero set (or the nodal set) of /. Then, SG \ Z{f) can be partitioned into 
finitely many connected domains Di,D2, ■ ■ ■ D^^, where Uk depends on the degree of 
the polynomial /. These domains {De}'^^^ are the nodal domains of /. We refer to 
[19] for color graphs that support the above bound on Uk. 



28 KASSO A. OKOUDJOU, ROBERT S. STRICHARTZ, AND ELIZABETH K. TULEY 

Remark 4.3. As mentioned in the Introduction, we have generated many data for all 
the OP constructed here. These data as well as the codes used to generate them are 
available at [19]. Certain data that are not included in the present article, deal with 
the dynamics of the OP at certain points of SG. That is, a choice of OP type and a 
fixed point x G SG yield a set of points {{pk{x) , pk+i{x))} in the plane. The passage 
from A; to /c + 1 is thought of as a dynamical system in the plane with some " attractor" 
toward which the points tend as /c — )■ oo. With our numerics, we had hoped to get 
a "snapshot" of this attractor. The computational complexity associated with our 
construction, limited us to to generate data for these dynamics only for 1 < A; < 200. 
Thus, we could not make some conclusive statements about the long term behavior 
of these dynamics. However, we refer the interested reader to [IS] for graphs of some 
of these dynamics. Our investigation in this context is parallel to recent results for 
the dynamics of certain OP related to self-similar measures on the unit interval [7]. 
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